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This is a short review of some RHIC results that have been most important for the small x physics community. 
We discuss saturation effects in deuteron-gold collisions, particle production in gold-gold collisions and some 
effects of the large "glasma" field configurations in the early stages of the collision. 



1. RHIC experiments 

The Relativistic Heavy Ion colhder (RHIC) at 
Brookhaven National Laboratory (BNL) (see e.g. 
Ref. [1] ) has been colliding gold and lighter nuclei 
at an energy of = 200A GeV since 2000. The 
accelerator also delivers polarized proton beams 
at a/s = 500 GeV, used to study spin physics. 
The main goal of the experimental program is the 
discovery and study of the Quark Gluon Plasma 
(see e.g. Refs. [5] for overall reviews of the exper- 
imental results). 

There have been four major RHIC experiments, 
the two smaller ones now decommissioned and the 
two larger ones being upgraded and taking more 
data: 

PHOBOS is a small detector, mostly of silicon. 
It has a huge acceptance in pseudorapidity rj 
and azimuthal angle 0, but mostly without 
Pt measurement or particle identification. 

BRAHMS is also a small detector, consisting of 
2 moveable spectrometer arms that give it a 
large rj coverage with particle identification, 
but a small acceptance. 

PHENIX is a large detector optimized more for 
electromagnetic signals (photons, leptons) 
with a good acceptance but not full az- 
imuthal coverage (A^ = tt in two spectrom- 
eter arms). 
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STAR resembles most closely a typical high en- 
ergy experiment, being built around a large 
TPC covering the full azimuthal angle. It 
is more optimized towards hadronic and 
global observables. 

The near term future of the BNL program is to 
continue taking data and improving the PHENIX 
and STAR detectors in parallel to the LHC heavy 
ion program. In the longer term there are plans 
to build an Electron-Ion Collider (EIC) to collide 
electrons with high energy protons and nuclei [3] . 
At CERN the LHC wiU, in addition to the pro- 
ton beams, have lead-lead (and eventually other 
species) collisions at ^/s = 5500A GeV. There 
is one detector (ALICE) dedicated to heavy ions 
and smaller heavy ion programs in the CMS and 
ATLAS experiments. 

An important tool for understanding the ge- 
ometry of the collision goes by the name of 
Glauber modeling. Experimental data are usu- 
ally presented divided into centrality bins deter- 
mined mostly by the total charged multiplicity. A 
Glauber model treats the nucleus as a collection 
of loosely bound nuclcons distributed according 
to the charge density distribution. These nucle- 
ons then scatter independently with the exper- 
imentally measured proton-proton cross section. 
A Glauber model the associates each impact pa- 
rameter with a certain number of participant nu- 
cleons A^part and binary collisions A^bin- One thus 
obtains an estimate for the typical impact param- 
eter corresponding to the centrality bins. While 
the assumption of independent collisions is not a 
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Figure 1. BRAHMS results on suppression of charged hadron production in the deuteron fragmentation 
region [4]. The plot shows Rcp (see text) for charged hadron spectra between two central bins (0-20% 
and 30-50%) and peripheral collisions (60-80% most central). 



very realistic one at high energies, Glauber mod- 
eling has turned out to be a convenient way to 
parametrize the geometrical aspects of the colli- 
sion; experimental data are usually expressed in 
terms of iVpart and iVbin- Experimental results for 
nuclear modifications to the spectra of produced 
particles are often presented as ratios (such as 
Raa, ratio of nucleus- nucleus to proton-proton; 
RdAu, deuteron-gold to proton-proton and Rcp, 
central to peripheral) corrected by the geomet- 
rical factor deduced from a Glauber calculation. 
These ratios are defined in such a way that in the 
case of independent scatterings of the nucleons in 
the nucleus they should be equal to one. 

We can understand the time evolution of the 
collision process in different stages: 

1. The initial condition at r = depends on 
the properties of the nuclear wavefunction 
at small x. 

2. Thermal and chemical equilibration. 

3. The Quark Gluon Plasma, surviving for 
some fermis around tq ^ t < 10 fm. If the 
system reaches local thermal equilibrium, fi- 
nite temperature field theory and relativis- 
tic hydrodynamics can be used to describe 
its behaviour. 

4. Finally, for r > 10 fm the system hadro- 
nises and decouples. 

In the remainder of this talk we will concentrate 



on the physics of the first stage, the initial condi- 
tion. 

2. Color glass and glasma 

The central rapidity region in high energy colli- 
sions originates from the interaction of the "slow" 
small X degrees of freedom, predominantly glu- 
ons, in the wavefunctions of the incoming hadrons 
or nuclei. At large energies these gluons form a 
dense system that is characterized by a transverse 
momentum scame Qs, the saturation scale. The 
degrees of freedom with pT ^ Qs are fully non- 
linear Yang-Mills fields with large field strength 
^ l/g and occupation numbers ^ l/cts, they 
can therefore be understood as classical fields ra- 
diated from the large x partons. Note that while 
this description is inherently nonperturbative, it 
is still based on weak coupling, because the clas- 
sical approximation requires l/ag to be large and 
therefore Qs ^ Aqcd- A "pocket formula" for 
estimating the energy and nuclear dependence of 
the saturation scale is Qs^ ~ A^/'^x~^''^: nonlin- 
ear high gluon density effects are enhanced by go- 
ing to small X and large nuclei. Ideally one would 
like to study the physics of the CGC at the Elec- 
tron Ion Collider [3], but already based on fits 
to HERA data and simple nuclear geometry we 
have a relatively good idea of the magnitude of 
Qs at RHIC energies [6]. The CGC is a system- 
atic effective theory (effective because the large 
X part of the wavefunction is integrated out) for- 
mulation of these degrees of freedom. The term 
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Figure 2. DGLAP parametrizations of nuclear pdf's from Ref. ^5j. The left panel shows the suppression 
of the nuclear gluon distribution required to simultaneously fit the nuclear DIS data from NMC (center 
panel) and charged hadron production data from BRAHMS (right panel) and PHENIX data. 



glasma [7] refers to the coherent, classical field 
configuration resulting from the collision of two 
such objects CGC. 



3. Deuteron-nucleus collisions 

A theoretically more controlled way to study 
the small- a; nuclear wavefunction is using proton- 
nucleus (in practice deuteron-nucleus) collisions, 
especially at forward rapidity, i.e. in the fragmen- 
tation region of the proton. In this case final state 
interactions in the plasma phase are expected to 
be absent, and the relevant degrees of freedom 
in the proton or deuteron are the relatively well 
understood large x parton distributions. 

The effects of parton saturation are particu- 
larly visible in the BRAHMS data (Fig. [1]) for 
charged hadron production, which shows a signif- 
icant depletion towards the deuteron fragmenta- 
tion region. Although these data DGLAP-based 
[5] analysis can reproduce this data, see Fig. [21 
this requires that the gluon distribution at small 
X is suppressed by a factor of ^^10 compared to 
the proton. This is a clear signal of the break- 
down of the DGLAP picture and the importance 
of saturation effects. The suppression pattern at 
large rapidities was predicted from the CGC [8] 
and later more detailed calculations have quite 
successfully reproduced it, see Fig. [31 
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Figure 3. CGC calculation of charged hadron 
production [9]. 



4. Particle production and geometry 

In order to compute particle production in the 
Glasma one is presented with the following sit- 
uation [1^. The valence- like degrees of freedom 
of the two nuclei are represented by two classical 
color currents that are, because of their large lon- 
gitudinal momenta (p^) well localized on the light 
cone (in the variables conjugate to p^, namely 
x^): ^ 5{x^). These then generate the clas- 
sical field that one wants to find. Working in 
light cone gauge the field in the region of space- 
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Figure 4. Classical field configurations. 



time causally connected to only one of the nuclei 
(regions (1) and (2) in Fig.|4]) is a transverse pure 
gauge, independently for each of the two nuclei. 
These pure gauge fields then give the initial condi- 
tion on the future light cone (r = \/2x+x~ = 0) 
for the nontrivial gauge field after the collision 
(region (3) in Fig. [4]). The field inside the future 
light cone can then be computed either numeri- 
cally [llj or analytically in different approxima- 
tions (see e.g. [T2] for recent work). The obtained 
result is then averaged over the configurations of 
the sources with a distribution W^j^[J''] that in- 
cludes the nonperturbative knowledge of the large 
X degrees of freedom. The resulting fields are then 
decomposed into Fourier modes to get the gluon 
spectrum. This is the method that we will refer 
to as Classical Yang-Mills (CYM) calculations. 

In the limit when either one or both of the color 
sources are dilute (the "pp" and "pA" cases), the 
CYM calculation can be done analytically and 
reduces to a factorized form in terms of a convo- 
lution of unintegrated parton distributions that 
can include saturation effects. Although this ap- 
proach (often known as "KLN" after the authors 
of [13] ) is not is not derived for the collision of two 
dense systems, it has advantage of offering more 
analytical insight and making it easier to incor- 
porate large-x ingredients into the calculation. 

The CYM calculations JJLj of gluon production 
paint a fairly consistent picture of gluon produc- 
tion at RHIC. The estimated value Qs ~ 1-2 GeV 
from HERA data [B] (corresponding to the MV 
model parameter g^fi « 2.1 GeV [TS]) leads to 
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Figure 5. Centrality dependence of the multi- 
plicity from a CYM calculation [14 . 



« 1100 gluons in the initial stage. Assum- 
ing a rapid thermalization and nearly ideal hy- 
drodynamical evolution this is consistent with the 
observed ~ 700 charged (~ 1100 total) particles 
produced in a unit of rapidity in central collisions. 

Of the more detailed geometrical aspects of 
the initial condition the basic features (such as 
the closeness to Vpart scaling) of the centrality 
dependence of particle multiplicities are mostly 
straightforward consequences of the proportion- 
ality of the multiplicity to Qs^] they are success- 
fully reproduced by both KLN and CYM calcu- 
lations [13114116] . see Fig.[51 

A striking signal of collective behavior of the 
matter produced at RHIC is elliptic flow; where 
the initial spatial anisotropy in a noncentral col- 
lision is transformed to an anisotropic distribu- 
tion of produced particles in momentum space 
through collective flow caused by the pressure 
gradients. Comparing hydrodynamical calcula- 
tions with measurements is a way to address 
fundamental properties of the medium, such the 
equation of state or viscosity, but this requires un- 
derstanding of the initial conditions. These initial 
conditions have traditionally been phenomenolog- 
ically parametrized in terms of the number of nu- 
cleons participating in the collision. The CGC 
provides a method to compute them from first 
principles, and has changed the interpretation of 
the experiments in terms of the viscosity [14117] . 
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Figure 6. Left: Two particle correlation measurement from STAR, showing the "ridge" structure elon- 
gated in pseudorapidity. Right: Long range correlations in multiplicity: b = ^ ^"^^^1 ''^^^^^^^2^^ ^ ; where F 
and B are pseudorapidity bins separated by A-q. 
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Figure 7. Parity violating correlation between 
reaction plane and momentum [19] . 



An example of a more exclusive probe of the 
medium is the J/^*. It has been called the "ther- 
mometer of the QGP" , because it the state is ex- 
pected from lattice calculations to melt at a tem- 
perature close to the QCD phase transition. The 
production mechanism of the J/^ from perturba- 
tive QCD is not very well understood, and since 
the charm quark mass is close to the typical val- 
ues Qs, it can be expected to be very sensitive to 
saturation physics [TS] , 



5. Glasma physics 

Correlations over large distances in rapidity 
can, by causality, only originate from the earli- 
est times in the collision process, and are thus 
sensitive to the properties of the glasma phase of 
the collision. Examples of such phenomena are 
the elongated "ridge" structure seen in two par- 
ticle correlations and long range correlations in 
multiplicity (see Fig. [S]) [20]. The boost invari- 
ant nature of the Glasma fields naturally leads to 
this kind of a correlation, and deviations from it 
should be calculable from the high energy evolu- 
tion governing the rapidity dependence Y21\ . 

Another remarkable phenomenon that is possi- 
ble in the Glasma field configuration is the gen- 
eration of a large Chern-Simons charge and thus 
parity violation |22| due to the nonperturbatively 
large field configurations. Through the so called 
"chiral magnetic effect" this can manifest itself 
in a parity violating correlation between the elec- 
tric dipole moment (or momentum anisotropy be- 
tween negative and positive charges; a vector) and 
the reaction plane (the positive charges of the ions 
generate a magnetic field perpendicular to the re- 
action plane; a pseudovector) . There are some 
preliminary indications in the data of a nonzero 
value for such an observable [T^j . 

In conclusion, experiments at RHIC have 
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clearly produced a strongly interacting decon- 
fined form of matter, but the full understanding of 
the experimental results requires an understand- 
ing of the small x physics giving the initial condi- 
tions of the collision process. This is the domain 
of strongly nonliear gluon fields at high energy, 
best understood in terms gluon saturation and 
the Color Glass Condensate effective theory. This 
framework provides a way to calculate properties 
of the initial condition from first principles, re- 
late them to DIS measurements, and understand 
some of the observables that directly probe the 
initial condition. With the much greater collision 
energy of the LHC these effects are likely to be- 
come even more prominent. 



REFERENCES 

1. J. L. Nagle and 
arXiv:nucl-ex/0203007, 



Ullrich, 



3. 



STAR, J. Adams et al, Nucl. Phys. 
A757, 102 (2005), arXiv:nucl-ex/0501009|; 
PHENIX, K. Adco x et al.,_ Nucl. Phys . 
A757, 184 (2005), [arXiv:nucl-ex/0410003l ; 
BRAHMS, I. Arsene et al, Nucl. Phys. 
A757, 1 (2005), arXiv:nucl-ex/ 0410020|; 

B. B. B ack et al, Nucl. Phy s. A757, 28 
(2005), arXiv:nucl-ex/0410022|. 

A. Deshpande, R. Milner, R. Venugopalan 
and W. Vogelsang, Ann. Rev. Nucl. Part. Sci. 
55, 165 (2005), |arXiv:hep-ph/0506148 . 
BRAHMS, I. Arsene et al, Phys. Rev. Lett. 
93, 242303 (2004), jajXiv:nucl-ex70403005l . 
K. J. Eskola, H. Paukkunen and 

C. A. Salgado, JHEP 07, 102 (2008), 
|arXiv:0802.0139 [hep-ph]l . 

H. Kowalski, T. Lappi and R. Venu- 
gopalan, Phys. R ev. Lett. 100, 022303 
(2008), (arXiv:0705.3047 [hep-ph]|. 



7. 



T. Lappi and L. McLerran, Nucl. Phys. 
A772, 200 (2006), |arXiv:hep-ph/0602189 . 
D. Kharzeev, Y. V. Kovchegov and 
K. Tuc hin, Phys. Rev. D68, 094013 
(2003), [arXiv:he p-ph/0307037 . 
A. Dumitru, A. Hayashigaki and J. Jalilian- 
Marian, Nucl. Phys. A770, 57 (2006), 



|arXiv:hep-p h/0512129|. 
10. A. " Kovner, L. D. 



H. Weigert, Phys. Rev. D52, 3809 (1995), 
|a rXiv:hep-ph/950532 . 

11. A. Krasnitz, Y. Nara and R. Venu- 
gopalan, Phys. Rev. Lett. 87, 192302 
(2001), I arXiv:hep-ph/010 8092| ; T. Lappi, 
Phys. Rev. C67, 054903 (2003), 
|arXiv:hep-ph/0303076'; A. Krasnitz, 
Y. Nara and R. Venugopalan, Nucl. Phys. 
A727, 427 (2003), |arXiv:hep-ph/0305112| . 

12. J. -P. Blaizot and Y. Mehtar-Tani, 
Nucl. Phys. A8018, 97 (2009), 
|arXiv :0806.1422 [hep-lAi] . 

13. D. Kharzeev and M. Nardi, Phys. Lett. 
B507, 121 (2001), | arXiv:nuc l-th/0012025|; 
D. Kharzeev and E. Levin, Phys. Lett. B523, 
79 (2001), farXiv:nucl-th/0108006|. 

14. T. Lappi and R. Venugopalan, Phys. Rev. 
C74, 054905 (2006), |arXiv:nucr-th/060902ll . 

15. T. Lappi, Eur. Phys. J. C55, 285 (2008), 
|arXiv:0711.3039 [hep-ph]|. 

16. H.-J. Drescher and Y.' N ara, Phy s. Rev. C 76, 
041903 (2007), [arXiv:0707.0249 [nucl-th]). 

17. T. Hirano, U. W. Heinz, D. Kharzeev, 
R. Lacey an d Y. Nara, P hys. Lett. B636, 

299 (2006), 

Drescher and Y. Nara, 



[arXiv:nucl-th/0511046|; 

Phys. Rev. 



H. J. 
C75, 



and 



034905 (2007), [arXivmucl-th/Oei lOlTl . 

18. D. Kharzeev, E. Levin, M. Nardi 
K. Tuchin, arXiv:0808.2954 [hep-ph], 

19. STAR, ^"S"^ A^ Voloshin, 



.arXiv:0806.0029 fnucl-ex 



20. J. Putschke, J. Phys. G34, S679 (2007), 



arXiv:nucl-ex /0701074 ; STAR, M. Daugher- 
ity, 'arXiv:0806.2121 [nucl-ex]; STAR, B. K. 
Srivastava, Int. J. Mod. Phys. E16, 3371 
(2008), a rXiv:nucl-ex/ 0702054|. 

21. A. Dumitru, F. Gelis, L. McLerran and 
R. Venugopalan, Nucl. Phys. A810, 
91 (2008), |arXiv:0804.3858 [ hep-ph]] ; 
F. Gelis, T. Lappi and R. Venugopalan, 
arXiv:0810.4829 [hep-ph] 

22. D. Kharzeev, A. Krasnitz and R. Venu- 
gopalan, Phys. Lett. B545, 298 
(2002), [ arXiv:h ep-ph/0109253 ; D. E. 
Kharzeev, L. D. McLerran and H. J. War- 
ringa, Nucl. Phys. A803, 227 (2008), 
[arXiv:0711.0950 [hep-ph] |. 



McLerran and 



